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AXISYMMETRIC SUPERSONIC FLOW IN ROTATING IMPELLERS! 


By Автнув W. GOLDSTEIN 


SUMMARY 


General equations are dereloped for wentropie, frictionless, 
axisymmetric flow in rotating impellers with blade thickness 
taken into account and with blade forces eliminated in favor of 
the blade-surface function. It às shown that the total energy of 
the gas relative to the rotating coordinate system 18 dependent 
on the stream function only, and that 4f the flow upsiream of 
the impeller is vortez-free, a velocity potential exists which is а 
function of only the radial and axial distances in the impeller. 

The characteristic equations for supersonic flow are devel- 
oped and used to investigate flows in several configurations in 
order to ascertain the effect of variations of the boundary condi- 
tions on the internal flow and the work input. Conditions 
varied are prerotation of the gas, blade turning rate, gas velocity 
at the blade tips, blade thickness, and sweep of the leading edge. 


INTRODUCTION 


No method of practical importance has been described in 
the literature for computing the flow in the three-dimen- 
sional region between the blades of a supersonic compressor. 
This deficiency in design technique has been partly met by 
using two-dimensional theory to compute the variations of 
the relative flow from blade to blade, by assuming no radial 
flow components, and by neglecting spanwise variations. 
This procedure is probably satisfactory for designs with а 
ratio of blade-root radius to blade-tip radius close to unity 
and for cylindrical hub and casing shapes. When this ratio 
is small, the spanwise variations are significant. A method 
devised at the NACA Lewis laboratory during 1950 is de- 
scribed herein for computing the magnitude of these varia- 
tions, when blade-to-blade variations are neglected, by 
assuming а large number of blades in the impeller and deal- 
ing with the circumferential average flow. This technique, 
when combined with estimates of blade-to-blade variations, 
is useful for designs with closely spaced blades. 

The usefulness of this approach may possibly be extended 
to supersonic flow in a manner analogous to the method 
described in reference 1 for subsonic flow; the axisymmetric 
solution is first found and then used as a base for computing 
the blade-to-blade variations. The particular method used 
in the blade-to-blade calculation of reference 1 is series 
expansion, which is not applicable to supersonic flows be- 
cause of the existence of solutions with discontinuous 
derivatives. 

In reference 2, the equations for internal flow with zero 


blade thickness are translated into the characteristic equa- 
tions for the derivatives of the stream function, and it is 
shown there that the equations are hyperbolic if the relative 
velocity is supersonic, provided the ratio of the tangential to 
the axial velocity is a given field function. (This condition 
includes the special class of blades containing radii from the 
axis of rotation.) If the tangential velocity is prescribed, 
the meridional velocity component is the decisive factor in 
determining whether the equations are hyperbolic. 

The present report treats the velocity components directly 
after the manner of reference 3 for arbitrary preassigned 
blade shapes. The design problem is solved by finding the 
hub shape that will give a prescribed velocity distribution 
along the casing (or any other selected stream surface). In 
the course of the solution, the tangentially averaged flow is 
completely determined in the interior of the impeller. 

The purpose of this paper, in addition to describing the - 
method in a mathematical outline, is to examine the per- 
formance of the technique in computed examples and to find 
whether any peculiarities in the flow arise from the boundary 
conditions imposed, and to develop a background of experi- 
ence for the design of impellers of the type considered by 
examination of the effect of systematic variations of impeller 
design parameters on the internal flow. The parameters 
varied are inlet prerotation, blade turning rate, gas velocity 
distribution at the casing, blade thickness, and sweep of the 
blade leading edge. 

Approximately 160 hours of computing are required to 
work out about 100 points for each of the impellers. This 
time can be reduced by application of the punched-card 
computing technique to part of the procedure. 


DEVELOPMENT OF EQUATIONS FOR AXISYMMETRIC FLOW 


Equations describing the fluid flow are developed in con- 
venient form by use of the equations of motion and con- 
tinuity, the conditions of axial symmetry, and isentropic 
state changes. The uniformly distributed blade forces, 
which are inconvenient in calculations, are eliminated from 
the equations by substituting the blade surface function, so 
that.the condition of integrability is automatically satisfied. 


EQUATION OF CONTINUITY 
The continuity equation for steady flow is 


V-(poW)—0 


t Supersedes NACA TN 2358, *Axisymmetria Supersonic Flow in Rotating Impellers" by Arthur W. Goldstein, 1951. 
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(Symbols are defined in the appendix.) If axial symmetry 
is assumed, the equivalent in cylindrical coordinates is 


2 (р orto (pwr)=0 


Tho presence of blades serves to modify these equations if 
the space filled by the blades is taken into account. If the 
impeller contains № blades with thickness ¢ measured in the 
direction of rotation, the ratio of volume available for fluid 
flow to the actual volume of. any space element approaches 
(n the limit as the number of blades increases indefinitely 
with Nt fixed) 
Á Nt 


I 


r 2ar 


The continuity equation is therefore obtained by reducing 

the specific-mass-flow terms ре and pw by the ratio A/r to 
account for the area reduction: 
д 

2. (04). (pwA)=0 (1) 


Because none of the scalar quantities is a function of @, 


equation (1) implies the existence of 8 stream function 


Vea such that ' E 
ф;=ршА 
¥.=—pvA (2) 
V1—0 
(subscripts indicate partial derivatives) or 
Vy=pAjx Y-—pAXW (2a) 


EQUATION OF MOTION 


Substitution of energy for pressure.—The equation of 
motion for steady, nonviscous, шеншоре flow in terms of the 
relative velocity is | ж э. Мо оё 


vL Wxvxwd тр 2WXO--F 


where wR is the centrifugal force and 217XQ, the Coriolis 
force. For isentropic conditions, the substitution 


1 
= uum Ь 
Р Vp-—V 
may be made, so that the equation is then 


v( MT - =W(vX W +20) +F=WXyXV+F 
(3) 
The — 


S hi arut) A 





is & measure of the total energy of the gas particles in the 
relative system of coordinates and remains constant (as will 


be shown for frictionless flow) for any particle; energy is 
interchanged among the forms of cnthalpy, kinetic energy, 
and potential energy (— w*r*/2) resulting from the centrifugal 
force field. In reference 2, Æ is shown to be constant for 
any particle for frictionless flow and no heat transfer with a 
finite number of blades in the impeller. 

In order іо show that Е is constant for any particle in 
axially symmetric flow, the direction of the field force is first 
identified. If the blades are considered to increase indefi- 
nitely in number and to retain their camber-line shape, the 
forces exerted by the blades arc merely the pressure forces 
exerted normal to the blade surfaces, which in the limit 
approach uniformly distributed forces in the body of the 
fluid. The relative flow lines are always contained in onc 
of the family of blade surfaces, so that the field or blado 
forces must be normal to the velocity. If the scalar product 
of equation (8) is taken with the relative velocity, the result 
is then 

W-vE-—O0 


which indicates that Æ does not change along any flow line. 
Elimination of blade force.—The blade force can be 
eliminated from the equation of motion if the equation for 
the blade surfaces is assumed known. The force may then 
be eveluated after the velocity components are known. 
Under the condition of axial symmetry, the equation for 
the family of blade surfaces 


a(r, 6, 2)=constant 


may be simplified. If the constant is changed, there is à 
shift from one blade to another. Because of axial symmetry, 
this shift is equivalent to rotating the original blade by a 
fixed angle. 
put into the form 

=ð +f (r 2) l (5) 


Because the relative velocity flow lines are contained in the 


blade surfaces, Ў із normal to Ve as well as to Vy. The 
relation among the three vectors can therefore be shown 
to be 


Wer V XVa (0) 


From the condition of axial symmetry, 


oH 
32 = =rjp-VE=0 

Because both VE and Vy are normal to j and ЇЇ, they 
must be parallel, and consequently Æ is known io be a 
function of y only (reference 4); that is, 


E-E() ! 


7 
vE=w S. ч) 


Therefore, the equation for the surfaces can be | 
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Insertion of equations (4), (6), and (7) mto equation (3) 
yields (without resolution of the term V X) 


= w= Vj. VX V) va- Va- YX v) Vy+EF 


The blade pressure forces F are normal to the blades and 
therefore parallel to Уе; this relation thus includes only two 
nonparallel vectors Va and VV, so that the components must 
separ&tely satisfy the relations 


F= (4 vj-VX v) Va—-—(rj.WxVvXV)va (8) 
and 
dE T $ 


Special conditions..—Two special conditions of technical 
importance are considered: In high-speed machines, blades 
with radial elements are desirable in order to eliminate bend- 
ing loads resulting from centrifugal forces. Another fre- 
quently occurring circumstance is that of zero vorticity 
upstream of the impeller. 

For blades with radial elements, the vector ¢ lies in the blade 
surface and therefore 

4-Va=f,=0 


Therefore f is a function of z alone. From equation (6), 


0— W.Va—- of twf, (10) 
so that for radial blade elements 
u=—wr 2i (10 а) 


This relation provides considerable simplification of the equa- 
tions developed later for solving the flow problem. 

For vortex-free flow upstream of the impeller, equation (3) 
is reduced for conditions there by setting Q—/'—VvxY-—0 
and substituting V for W to get 


©(%+-5-)=0 





or 


h-+~=constant 


Also, zero vorticity gives, for axial symmetry, 


0=VXV=—1 = (utor)t7 2: See 2 (ur- ет?) 


2 т Or 
From the first and last terms and axial symmetry, 


ur--uwr*-constant 


Then the value of the total energy F relative to the impeller is 
72 2422 үз 


which is constant upstream of the impeller. Then by equa- 
tion (7), E remains constant in the impeller,so that the equa- 
tion of motion (9) becomes 


Va-VXV=0 (11) 
If the system of coordinates used is г, z, and æ for expressing 
equation (11), there results (reference 4) 


E (V-vexk)| =| > (rV. ixva] 


The subscripts indicate which variables are held constant. 
This equation and the fact that the quantities in parentheses 
are functions of r and z only establish the existence of а poten- 
tial function of 7 and z only ($(r,z)) such that 


E —rVy-.Voxk 
or 

Ф -— 
FI XVa 
oe] _ 

де, 


The last relation indicates that with r and z constant, Ф is 
constant; this is true regardless of the third independent 
variable, which may be taken as 0. The three equations for 
$ are then put into vector form 


vé—i(rV-Vax E) J-k(rV-iX Va) (12) 
'The absolute velocity may then be shown to be 
V=VO-+-r(utor\Va (13) 
and from 
F -Va=W-Vatourj -Va=o 


equation (13) gives for the moment of momentum y, 


р erben (14) 
so that 
=vep oe Va (13a) 


The potential Ф can be used as an alternative to the stream _ 


function in solving for the flow. Equation (13) indicates 
that the projection of УФ on the blade surface is equal to 
the projection of the absolute velocity on the blade surface. 
The equation for y and that for Ф are obtained in the following 
manner: The velocity V is eliminated from equation (9) by 
means of equation (6) and there is obtained 
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Va: vx(7 VUXVa)+2ufe=—P or 


which is also applicable in the general case with nonzero 
vorticity upstream of the impeller. 

In order to obtain an equation for d, equation (13a) is 
substituted into equation (2a) 


Vi — pArVOX V=pArvexVe-t pAprvéX Va 
Then 
0— Vx V) VX [pArv8X (УФ uVo)] 


Although obtaining a solution in terms of y is feasible 
whether dE/dy is zero or not, when dE/dj —0, Ф may prove 
more convenient. In regions where W'?/aà? is nearly 1.0, 
because |W| is а double-valued function of the stream- 
function derivative e| W |, it will be uncertain whether the flow 
is subsonic or supersonic. А solution obtained from а net- 


work would require а, very fine spacing in regions where 


DT | ea. 
GENERAL EQUATIONS IN CYLINDRICAL COORDINATES 


In general, where f, 740 and dE/dy #0, the equations for the 
internal flow are solved simultaneously for the two velocity 
components v and w, when u has been eliminated from the 
equations by using equation (10). The equations are trans- 
lated from their vector form to scalar notation by utilizing 
cylindrical coordinates. Then а further transformation is 


employed with characteristic quantities: as the а 


variables. 


The continuity equation (1) is first written in vector form. 


(the condition of axial symmetry being used) and multiplied 
by r/pA 


IU G 4 w)=v. Wy: EE W -V(A/r)=0 


For isentropic flow, 


Introduction of E by equation (4) and use of W’-VE=0 
produce 


W. унде „дА ү 


ти 
Е а? а? +30 E S e : 


'This expression finally reduces to 


(1 yu Ho) (175) 309) 
(п) ORA wi lg ^ qs 


In cylindrical coordinates, equation (9) becomes 
dE 
0, — ved fts fa pA | ы (16) 


If equation (10) is now substituted into equations (15) and 
(16), there resulte for continuity 
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[m EA), o (“= —® Eum w(t- a2) ot 


(п) o Ef fap 200 faf 

: 2 
"Gem e aante anjo ат 
and for the equation of motion 


(1 -r5f2)w, +r fL f (0, —w)—(1 I 


njal o rez -2 f AT (18) 


USE OF CHARACTERISTIC VARIABLES 


Statement of basic equations.— The method of reference 3 
is applied to equations (17) and (18). "The objeet is to find 
two alternative equations with new variables and », instead 
of z and r, such that each equation (corresponding to equa- 
tions (17) snd (18)) contains derivatives with respect to only 
one of the independent variables. Two additional equations 
are required to find the position coordinates 2 and 7 in terms 
of the characteristic coordinates & and т. There are two 
characteristic slopes dr/dz, t4 and ¢-, the first corresponding 
to £ varying and у constant, the second to £ constant and з 
varying. If t represents either of these slopes and с repre- 
sents either of the coordinates corresponding to the appro- 
priate ¢, the results are 


Te™ {Ze (19) 


[9 уара та) e» 


(the upper sign of + refers to ¢,, the lower, to t-) and 


Jw, - G p, -- L2,—0 (21) 


where 
B—1 rfe yap 
М? (u? + v? Fw?) 


a? 
Ju рау 
G (M rf, f ere) EV NFA 
L=— (r ffit tl +r fA 
Tru f, —1)— (1-5 wt Ay |р, 
D, (wtf 2wv fet "у + 


"(1 ou) 49r 204/1) l/r) | wr o(Ajr) 
Á or TX Oz 


D,=rf Ades) ez) —2(u- er)fi— PASS 


(21а) 
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Equation (20) for ¢ shows that the system is hyperbolic if 
l. 

In the homogeneous equations (D, —D,—0), the coefficients 
Г. and L. (corresponding to ft, and f.) are zero; conse- 
quently, in the hodograph plane the slope of the characteris- 


tics would be given by $2=—% In plane flow, J/G is а 


function of v and w only, and this equation may be integrated 
and the hodograph solutions tabulated in advance of applica- 
tion to any particular problem that may arise. In the pres- 
ent circumstances, however, J/G also involves the position 
coordinates, and D,~0, 1550. A correspondingly easy 
method of solution by simple waves is therefore unavailable. 





= 
Fiavke 1.—Characteristic and cylindrical coordinates. 


Instead, the four (c—£ or 7) equations (19) and (21) are inte- 
grated simultaneously in a step-by-step procedure. 

Method of computation.—In figure 1, conditions are sup- 
posedly known at P; and P, and conditions are to be computed 
at Ps. Then equation (19) in difference form is 


Tg—Ti— (23— 21) f+ (22) © 


7з——Те== (Z3— 20) f— (22а) 
This system may бе solved for r, and 2. Then similarly for 
equation (21), the difference form is 


Jiws — w) + GL (05 — 91) + Ly (23— а) —0 (23) 
J (w — w) + G.. (v4 — v9) Е. (23 — 22) —0 (23a) 


The subscript , indicates that т is constant and the subscript 
— indicates that £ is constant. In equations (22) and (23) 
the coefficients t, J, G, L are averages between values for the 
points P,and P4. For the first approximation, only the values 
at Р; аге known and these values are used. 

When approximate velocity components and position coor- 
dinates are thus determined for Ps, the coefficients t, J, G, L 
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may be computed there and averaged for a better approxima- 
tion of the position coordinates and the velocity components 
at P, This process is repeated as often as necessary to 
arrive at an answer of satisfactory accuracy. The smaller 
the interval, the quicker the convergence to the solution. 
In each stage of the process, if dZ/dy/ 750, then ұз must be 
computed from 


Ае Í. РЧ) f КОЖЕ УГУ 


From the value of уз, 05/20 at Ру may be found. 





Ё constant 
————— 2] consfont 


FIGURE 2.—Relation between shape of bouudary and boundary conditions. 


BOUNDARY CONDITIONS 
GENERAL 


The discussion in the previous section implies that boundary 
conditions are known on & curve C containing points such as 
Р, and Р, and that the solution is then extended into а region 
containing points such as P, The shape of the blades and 
Е (y) are assumed known in that region. On the curve C, 
there are known in addition the values for w aud т, from 
which all other variables (such as 68/40, v, u, p, a?, and so 
forth) may be computed. Suppose the initiel curve is such 
that the characteristic curves intersect it at two points, such 
as А and D ог A and E in figure 2. Then by equation (21), 
there is a relation between the variables at A and D that does. 
not permit the independent specification of the variables 
v and w at А and D. Similarly if two characteristics passing 
through two boundary points reach & third boundary point, 
then another degree of freedom is lost. In figure 2, if two 
variables are prescribed on AB and EF, then one variable may 
be prescribed on BC and DE, and none on CD. 


ТІР CONDITIONS 


The supersonic impeller has a high pressure ratio per stage 
and a large air-handling capacity per unit of frontal area. 
In order to maintain the large capacity, the tip radius at 
the entrance should be at least as large as at the exit. In 
order to provide a maximum work output with a minimum 
turning, the exit tip radius should be at least as large as the 
entrence tip radius. А reasonable shape for the tip stream- 
line is therefore & surface of constant radius, although smaller 
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radii between the entrance and the exit might be useful in 
order to affect the internal flow. 
lected, then, as it is also а streamline, the ratio v/w is 
specified. Selection of the velocity or of a velocity component 
remains. This selection may be made by consideration of 
the blade-to-blade velocity variations and by using some 
approximate method such as that outlined in the following 
paragraphs. | mS 

In three-dimensional motion (finite number of blades), 
Kelvin's theorem states that if the absolute velocity is irro- 
tational upstream of the impeller, it will remain so inside the 
impeller; that is, | 


VXV'zVxW'--2ek-—0 


where the prime indicates velocities of the nonsymmetric 
flow. Then 


Vy-Vx 770 
may be reduced to 
д T; = ou! f EC 
36 CW i-r Fe 1 T2wr)g4 


where Ow'/óe is the arc-length partial derivative of u’ 
along the direction of 9. This equation is now integrated 
from the pressure surface of one blade to the suction surface 


of the next in the direction r=constant, z—constant. The 


variable 6 increases by the amount 27A/Nr, so that - 


vj ук "8s Ou’ _.{%_, 2тА„ =. 
СИ) (И. p. 3: 2+7) u'd04- N 2079-1 





If it is now assumed that the values of 00/05 and ої u’ 
averaged across the channel are equal to the values obtamed 


A C 





A 
(а) 
NU _ | 
A Blades a 7 
А 
2 (b) 


fa) Conditions on B'B and BC unrelated. 
(b) Conditlons on B’B related to those on ВС, 
FiGURE 3.— Boundary conditigns for impeller. 


If the casing shape is se- - 


(24), 


= m» P ЕЕ ГУ TRAC : 
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from the axisymmetric solution, and further that (H"-¢),= 
IW], then 


| yr | т (52 a) та тА 


UNS (24u) 
if the relation 
of |=29 


is utilized. In this equation all the velocity components 
without subscripts are obtained from the axisymmetric 


solution. If 
ИЕ I—2] T] 


or some other connection between |0,1, (Ml, [JP], and 
the impeller geometry is assumed, the blade surface veloci- 
ties at the casing may be obtained with that relation and 
equation (24a). The distribution of |W'| at the casing may 
therefore be selected on the basis of obtaining a desired 
blade surface velocity distribution or distribution of v. 

Suppose that, in some manner, all the variables have 
been selected at the casing (curve BC in fig. 3(8). The 
flow is then uniquely determined in the region BP,C enclosed 
by the casing and two characteristic lines. If the flow is 
also prescribed at the leading-edge surface BB’, then the 
flow is determined in the extended region B'BCP,Ps. When 
the characteristic line BP, falls inside this region (fig. 3(8)), 
no difficulty results. If, however, the characteristic B'P, 
intersects the casing inside the blade region (fig. 3(b)), then 
the velocity may not be prescribed on the curve BP, (refer 
to discussion of fig. 2), or alternatively, only one of tlie 
velocity components may be prescribed on BB’ (BB? is not 
a streamline). Similar considerations are valid on the curve 
CC’ at the exit. 

The selection of values for the variables » and w on the 
entrance edge BB’ must be made by joining the solution 
for the flow in the region ABB’A’ with that for the flow in 
the blade region BCC’B’. Because the blades may deflect 
the air at BB’ discontinuously and because the effect of 
thickness may be felt by the air in a discontinuous manner, 
the discontinuities.in the flow should be evaluated. A 
method for making this estimate follows: 

The angle 8 is defined as the angle between the meridional 
projection of the blade leading edge and the z-axis (fig. 
3(а)). The conditions outside the blade region are desig- 
nated by subscript zero, those inside the blade region, by 
no subscript. Then across the surface BB’ the equation 
of continuity yields 


pos sin B— vo cos 9-(1 ~ >) p(w sin B—* cos В) (25) 


The discontinuity in blade thickness can cause no force 
in the entrance surface BB’, and the blades can exert no 
force in the blade surface. Therefore no force and no 
discontinuity in momentum exist in the direction ‘of the 
curves representing the intersection of these two surfaces; 
that is, 
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(W — W)-(t cos 8—£ sin 8)XVa=0 


or 
r(u—ug)(f, sin B4- f. cos 8)—(»— vo) sin B-+(w—wy)cos B (26) 
There is also available 

uo —r(w fdv fr) 


If the process is assumed isentropic, then 


On a 


is obtamed by utilizing the fact that 





2 252 = 243 


and that 
Tyg—r 

There are then available equations (25), (26), (27), and 
(28) for finding р, u, v, and w. The initial values of f, and 7, 
may be chosen from these equations to obtain desired rela- 
tons between the velocity components outside and inside the 
blade regions. 

APPLICATIONS OF METHOD OF CHARACTERISTICS 

In the examples computed in this section, the conditions 
for which solutions are obtained vary. The purpose of these 
solutions is to find the effect of some variations of prescribed 
conditions on the resultant flow. For all the examples, 
dEH/dt=0; the flow is vortex-free upstream of the impeller; 
at the leading and trailing edges the radial component » of the 
velocity is zero; and the blades contain radial elements 
(3.==0). The relative Mach number at the conjunction of 
the leading edge of the blade and the impeller casing is 1.94. 
Except in example VI, the value of the blade angle with 
respect to the axial direction is zero at the exit (f,—0). The 
other conditions are summarized in the following table: 








































Entrance 
| 
: 'Sweep-| Blade ; 
ing; Turning : m 
p Example | prenota. | Axial | ' rate | ba thick: 
І tion, fnlet 
| (u--er)r | velocity 
I 0 0.85 а, | 
| | 
j II —0.3 аш: | 0. 835 а, 
i 
| Ш —0. 3 аът: | 0. 835 a1 
IV ~—0. 3 arr, | 0. 835 а 
| —0. 3 а: | 0. 835 ay 
2250] 


у 
а  ——— —— 
VI | ~O.3 ar: 0. 835 a; | 1. 





Ы pend. by conditions at leading edge to z/r;—0.350; held constant at 2.1 a from 
че о 
b Constant at 1.94 ato z/ri» 0,25; decreased to L4 а (1.6 a4) at ext 

€ Zero from z:x0 to zjrc20.25; from z[T4-x0.25 to z/r, 0.5 the usi m for blade thickness is 


Amr exp iT [ 1-с oF (4-5) |} 


4 Determined by leading edge initia! value, LM a. 
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Example I.—The results of the computation for the first 
example are summarized in figure 4 (a), which shows stream- 
lines, contours of constant Mach number and constent energy 
level, and some of the characteristics of the £-family, which 
intersect the leading edge. Because of the intersection of 
the characteristics from the leading edge with the tips of 
the blades, there is no control over the Mach number, as the 
blade shape was prescribed by 


5097 df 1 .04 — Z (14-cos r=) 


The Mach number at the casing rises to a maximum of 2.1 at 
z[r,— 0.356, and is then assigned а constant value. By equa- 
tion (4), the stagnation enthalpy Н (2h4-V?/2), which 
measures the energy of the gas, increases with z at the rate 


oH Ou 


dz de 
At the casing v=0, and therefore, from equation (10а) and 
M= tHe, 


= cr 
WI 


rw? d?f 


ul Wl. TFT dz? 


This expression shows that the rate of work input is increased 
by increasing the rate of blade turning and decreased by 
increasing IF. Figure 4 indicates that the rate of turning of 
the blade is too small in the beginning of the impeller and 
that as a consequence nearly all the work input is concen- 
trated in the back portion of the impeller. This condition 
is ageravated by the acceleration of the gas in the initial 
portion (u[}¥}.<0 because u<(0). Increasing the turning 
rate, however, wil open the flow area and cause |I¥| to 
increase unless the area is controlled by blade thickness. 


A great disparity of work exists between the root and the 
tip; the values for AH fa? are 0.65 and 2.5 at the exit. For 
isentropic compression, the corresponding ratios of outlet to 
inlet stagnation pressures are 2.48 and 9.0 for the root and 
the tip, respectively. The Mach number at the root at the 
exit (z—r,, 7=0.507 rj) is in the neighborhood of 2.25, and 
а maximum value of about 2.5 is reached at z—7;, r—0.70 r;. 
Subsonic Mach numbers are reached in the region 
0<2/",<0.3. The solution could not be extended into this 
region. 

If the velocity were caused to decrease at the exit (2=7;,) 
by decreasing the velocity at the casing (r—7;) in the rear 
part of the impeller (0.356<2/r,<(1.0)}, the mass-flow density 
at the exit would increase with the decrease in velocity and 
the required flow area would decrease. Termination of the 
hub would be expected, not at 7—0.507 r, but at some larger 
value where the work output would be higher. Another 
method of reducing the blade height at the exit and thereby 
reducing the work variation there would consist in relaxing 
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(8) Example I. 
Figure 4.—Flow characteristics. 
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(b) Example II. 
Figure 4.—Continued. Flow characteristics. 
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the condition of 27—0 and utilizing some other relations that 
would correspond to lower Mach numbers. Two variables 
or relations could also be prescribed at the exit, but conditions 
on the blade tip (r—7;) could not be specified for some region 
near the back of the impeller. In this respect the entrance 
and exit conditions would be similar. 

Example II: Prerotation of gas.——In order to maintain 
supersonic flows throughout the impellers so that the flow 
could be computed, the relative entrance Mach number at 
the blade root for example TI was increased by assuming an 
initial counterrotation. The rotative speed had to be reduced 
from 1.68 a, to 1.44 а to maintain the work output, the 
relative Mach number at the leading-edge tip, &pproximately 
the same blade angle there, and the same axial velocity com- 
ponent at the inlet. This reduction had & further beneficial 
result of reducing the absolute Mach number at the exit for 
a given relative velocity axially directed. The condition of 
z—constant for the leading edge cannot be maintainer! for 
shockless entry with radial blade elements and a counter- 
rotating free vortex entering the blade system. The sweep 
angle 8 chosen for this example was 45°, enough to permit 
free selection of W—1.94 a over the entire blade tip, because 
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(c) Example ШП. 
FIGCRE £—Continued. Flow characteristics. 
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the characteristic curves from the leading edge do not inter- 
fere with boundary conditions on the casing. The reduction 
in velocity at the casing reduces the blade height at the exit 
and the work-output variation there. 

In order to examine deviations from an exemple more 
nearly representative of a good design, the velocity at the 
blade tip for example Н was not maintained the same as in 
example І, and consequently the differences between examples 
[ and П are not ascribable solely to inlet prerotation. Com- 
parison of the effect on rotative speed (for the same turning 
at the tip and the same work output) and on diffusion at the 
blade roots is, however, valid. 

The results of the computations are shown in figure 4 (b). 
In this example, the turning rate is predetermined in the 
region —0.5<2/r:<c0 by the velocity distribution upstream 
of the leading edge, the condition of radial blade elements, 
and the slope of the leading edge. The turning rate is con- 
tinuous at z— 0, and satisfies the equation 


£ s )=3.008—2.298 (32—1)} (0< z/ri« 1/3) 


—3.008[1—1/4(32—1)]  (1/3<2/r;,<1) 


532 


Compared with the first example, the work-input rate at the 
- blade tip is greater in the initial] section; however, the work 
distribution is still unsatisfactory. The exit blade height is 
also much less, with the result that the work input is much 
more nearly equalized; at the root the value of A/a? is 1.8, 
whereas at the tip the value is 2.5. These values correspond 
to stagnation pressure ratios of 5.2 at the root and 8.35 at the 
tip. This variation could be reduced by sweeping back the 
trailing edge to provide space for additional work mput at 
the hub streamline. This sweepback is more feasible with 
the smaller blade height and work variation of the second 
example than for the first example. 

Example III: Effect of turning rate.—In order to reduce the 
size of the impeller by increasing the turning rate, example 
III was computed with twice the blade turning rate of exam- 
ple II for 2720 and with the same over-all work output. One 
result is the greater diffusion on the hub, where a minimum 
Mach number of 1.2 is reached instead of 1.4 as in the pre- 
vious example (see fig. 4 (c)) The imereased diffusion 
probably results from the increase in hub curvature, which 
in turn is caused by the shorter axial distances required for 
the changes. Another result is the increased. height of the 
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blade at the exit, and an increase in the exit Mach number 
from 2.3 to 2.7. The increased streamline curvature also 


. accounts for the higher exit velocity and this higher velocity 


in turn explains the larger blade height. The characteristic 
line through the casing at the exit (r—7;, 2—0.5 rj) intersects 
the hub at 2=0.155 r, The characteristics through the re- 
maining portion of the hub (22>0.155 r;) are affected by con- 
ditions at the exit (2=0.5 r,) where the velocities increase 
with decreasing radius. The result is a large increase in 


velocity at the blade hub in the region 2>0.155 ra 


Example IV: Effect of diffusion at casing.—In the next 
example, all conditions remained the same as in the third 
example except that the prescribed velocity at the casing was 
varied smoothly from the value of 1.94 а at 2—0.25 r, to 
1.40 a &t.2—0.5 г. The desired reduction has been obtained 


‚ in blade height and work-output variation, but the discharge 


velocity has а large value of 2.6 a (fig. 4(d)). It may tbere- 


fore be concluded that the alternative of directly prescribing 
а condition with low relative velocities at the exit would have 
been more effective than prescribing conditions at the blade 
tip for reducing the blade height at the exit. 
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FIGURE 4.—Continued. Flow characteristics. 
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А new ph:22?menon has also arisen in that the character- 
istics passing through 2-0.5r, 7—0.92r, and 2-0.5r, 
r--0.808r, intersect at 2--0.397r, 7—0.832r, At this 
point, the Mach number changes discontinuously from 2.13 
to 2.43 and the streamline inclination to the axis changes 
from 14.9? to 7.4?. This single point.does not represent the 
full extent of this difficulty. Actually, there is a region 
bounded by & curve in which two solutions are obtained; 
each solution is determined by a different region of the 
boundary. Consequently, the boundary conditions are 
incompatible with obtaining a single-valued continuous 
solution. This situation is analogous to the Prandtl-Meyer 
flow, where the characteristics intersect at the corner. If a 
computation were started for a simple two-dimensional 
supersonic flow from a streamline shape and from the velocity 
distributed thereon, which were exactly the same as on one of 
the streamlines of a Prandtl-Meyer flow around a corner, 
then, when the solution was extended toward the corner, the 
characteristics would intersect at the corner in the manner 
of figure 4 (d). Тһе intersection on figure 4 (d) can be taken 
to represent a physical boundary to the fluid in that the 
streamline passing through that point would be taken for 
the hub shape. Only 90 percent of the flow (ф=0.9) could 
be accommodated in the impeller of example IV if continuous 
flow with the imposed boundary conditions were required. 
An obvious method of increasing the flow is to increase the 
slope of the characteristic lines at the exit for r<c0.9r;. This 
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increase could be accomplished by decreasing the relative 
velocities there and discarding the condition v—0 to obtain 
а shorter blade height. 

Example V: Increasing blade thickness.—In none of the 
examples computed so far was blade thickness considered. 
An obvious effect of increasing blade thickness is to decrease 
the velocity in a region where there is no independent con- 
trol of the velocity (such as the region 2< 0.357, in example I, 
fig. 4 (a)). This phenomenon provides the degree of freedom 
required in such regions to control the relative velocity || 
as well as the surface velocities |Р, and [W,|. In regions 
where the velocity |W | at the blade tips is under control (the 
entire tip sections for examples H, IIT, IV), altogether dif- 
ferent effects may be expected. Тһе effect of thickness 
under these circumstances was investigated by computing 
another example with a cosine-like bump in the blade thick- 
ness starting at g=0.25r, (See the previously presented 
table of boundary conditions.) The maximum thickness 
was reached at 2=0.375r; At r—r, the thickness is 20 
percent of the blade spacing, and at the hub the value is 25 
percent. The results of the computation are shown in figure 
4 (е). The effect of the bump is not felt upstream of the 
characteristic with positive slope passing through z—0.257,, 
r==r, The most obvious effect is a decrease of flow Mach 
number in the region of increasing thickness (and inward 
from the casmg, where the velocity is not preassigned), as 
might be expected from one-dimensional theory. At z= 





(е) Example V. 
FicugRE 4.—Continued. Flow characteristics. 
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0.25r, and 22 0.375r,, thechangeinbladeheight is about 10 per- 
cent аз compared with example IV. Another effect of 


thickness is to be seen in the waviness of the hub streamline | 


in the region 0.15<2/r,<(0.25 where the thickness has not 
yet started to increase. This bump may be traced in other 
streamlines along the characteristic running through the 
point of maximum thickness at the casing (r—r;, 2—0.375 т). 
The effect of the bump on the hub is therefore to increase the 
local velocity by curvature in order to hold the required 
velocity at the casing; this action from hub to casing takes 
place along the characteristic line as simple waves in two 
dimensions, and the velocity is seen to have increased all 
along this characteristic line. . The intersection of the char- 
acteristics and the double-valued solution occur near the 
exit over a considerably larger portion of the flow than in 
the previous example. In this example also, if a lower ve- 
locity were prescribed at the exit, a unique solution would 
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be obtained. It would also be feasible to carry out this 

idea and at the same time to prescribe another condition 

(such as »=0) provided that one were willing to lose control 

of the condition at the back portion of the casing (r=r) in 
a manner analogous to that used at the front portion in 

example I. 

Example VI: Effect of leading-edge sweep.—1t is of more 
technical interest, however, to carry out this scheme at the 
entrance, for then the leading-edge sweep may be reduced 
from 50 percent of the axial depth of the machine to consider- 
ably less; А computation of the inlet flow conditions for 
such an impeller is shown in figure 4 (f), where the conditions 
at the lesding edge are the same as for example III. The 
blade shape is different, with the second derivative of angle 
d'eldzi— —d'f(z)|d£! constant. Also, the velocity at r=1.0 
cannot be preseribed. Only in the region enclosed by the 
casing and leading-edge characteristics was the computation 
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FIGURE 4—Coneluded. Flow characteristics. 
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completed. The result (fig. 4 (f)) shows a large increase in 
Mach number, especially at the blade tip. Because of this 
increase, the blade turning accomplishes no work in the 
beginning of the impeller, but shows a region of energy 
abstraction from the air; the initial negative rotational com- 
ponent becomes more negative as the air turns axtally and 
the velocity increases as in example I. The situation is 
quite different in the back portion where the blade angles 
approach zero and the tangential component does likewise; 
all the work input is concentrated in the back portion of this 
impeller as in the other examples. . Postponement of work 
input to the back portion of the impeller is to be expected in 
any case where the velocity is increasing as the blades turn 
toward the axial direction, or where the velocity is decreasing 
as the bledes turn past the axial direction. Internal flow 
and work distribution could be considerably improved in 
this example by increasing blade thickness in the region of 
large velocity increase. 


CONCLUSIONS 


General equations developed for isentropic, axisymmetric 
fiows in rotating impellers with any given blade shapes 
showed that the blade-force term can be eliminated and 
solved for, if desired, after the flow solution is completed. 
Further, if the flow upstream of the impeller is vortex-free, 
then a potential function of only radial and axial distances 
exists, which may be used to find the projection of the 
absolute velocity on the blade surface. 

For supersonic flow, the method of characteristies was 
used to compute the flow for several examples with radial 
blade elements and vortex-free inflow. These examples 
show that 

(1) With zero radial inflow, the influence of the leading 
edge and the blade turning on the casing velocity in regions 
where this velocity cannot be prescribed in advance causes 
an increase in velocity and a reduction of blade loading and 
work input in the initial portion of the compressor. These 
effects can be compensated for to some extent by increasing 
blade thickness. 

(2) Sufficient leading-edge sweepback adds a degree of 
freedom to the boundary conditions at the blade tips and 


therefore permits a more constant or a decreasing velocity 
Inside the impeller. 

(3) Reduction in the exit velocity results in a shorter blade 
height and a more uniform distribution of work input at the 
exit. 

(4) Increase in the rate of turning causes more rapid hub 
curvatures and consequently larger internal velocity gra- 
dients. Diffusion on the hub is therefore greater in the 
initial part of the impeller, and a lower velocity minimum is 
reached there, whereas the gas velocity at the hub near the 
exit reaches larger values. This situation requires larger 
blade heights at the exit with a larger exit velocity at the 
root and greater variation of work output from root to tip. 
Diffusion of the gas velocity at the tip will moderate this 
effect, as will also a relaxation of the condition of zero radial 
velocity at the exit and the substitution of a condition giving 
lower exit velocities. 

(5) Not every set of boundary conditions results in a 
continuous, single-valued solution in the region of interest. 
In the examples considered, a relaxation of the condition of 
zero radial velocity at the exit would have permitted another 
condition giving lower velocities there. This condition 
would have eliminated the region of the double-valued 
solution and given a shorter blade and a more nearly 
equalized work output. 

(6) The effect of variations in blade thickness in regions 
where velocity is also prescribed at the blade tip is to change 
the velocity as in one-dimensional flow but not by the same 
amount, for there is some compensation by variations in 
blade height. The prescribed velocity and thickness require 
hub curvatures that cause velocity changes to be projected 
like simple waves on characteristic curves. It is therefore 


· possible to find regions of increased velocity when the blade 


thickness is increased. 


Lrwis Ётлент PROPULSION LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
CLEVELAND, Onto, January 24, 1951. 


APPENDIX 
SYMBOLS 


The following symbols ere used in this report: 
А т—А 2т 


а sonic velocity of gas 

Qi sonic velocity of gas at conjunction of leading edge 
and impeller case 

a* critical speed of gas upstream of impeller 

B an abbreviation (see equations (21a)) 


abbreviations (see equations (21а)) 


E sum of stagnation enthalpy and potential energy 
of gas (specific) in moving coordinates 


WE ir 
(=®+т—®) 
distributed blade force per unit mass of fluid (а 
vector) 
f part of blade surface function dependent on r and z 


only (a—8--f(r,z)) 


an abbreviation (see equations (21a)) 

absolute stagnation enthalpy (h-4-V?/2) 

work input. to gas (rise in absolute stagnation 
enthalpy) 

static enthalpy of gas (specific) 

unit vector parallel to curve. é=constant, z= 
constant (4=Vr) 


Rue g 





J an abbreviation (see equations (218)) 

j unit vector parallel to curve z—constant, r—constant 
(j—7V0) 

k unit vector parallel to curve r- constant, 0— 
constant (k=Vz) 

L ап abbreviation (see equations (21a)) 

РА 

М local relative Mach number ( T) 

N number of blades 

p static pressure of gas 

q unit vector parallel to И” 

R vector of 7 

r normal distance from axis of rotation 

"i distance from axis of rotation to blade tip at leading 
edge 

ds arc-length element parallel to W 

É blade thickness measured in directjon of rotation 

u rotational component of relative gs velocity 

y absolute gas velocity (a vector) [io+} (u+ ar)+ kw] 


radial component of gas velocity . 
W relative gas velocity (a vector) (o+jutkw) 
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w axial component of gas velocity 

2 distance measured parallel to axis of rotation | 

a function for blade surface (a (r, 0, z)=constant) 

8 angle between meridional projection of leading edge 
and axis of rotation | 


ү ratio of specific heats 

t either £4 or t.. in general form of equations 

t4 slope of curve for £ variable, y constant 

С: slope of curve for у variable, £ constant 

7 one of characteristic coordinates for gas flow in 
impeller 

0 meridional angle (for cylindrical coordinates) | 

u absolute moment of momentum (r (u+ wr}) 

t ` опе of characteristic coordinates for gas flow in.. 
impeller 

p . gssdensity = 2  1— 

т either 7 or £ in general form of equations 

Ф velocity-potential function of r апа z only that exists 


for initially irrotational absolute motion 
y stream function (VypÀj XV) 
Q vector form of w 
W angular velocity of impeller 


Subscripts: 

0 upstream of impeller 

p pressure surface of blade at tip 
8 suction surface of blade at tip 


Subscripts r, 6, 2, с, £, 7 indicate partial differentiation. 
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